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Abstract 

We show that the recently developed Hamiltonian theory for high energy evolution in QCD 
in the dilute regime and in the presence of Bremsstrahlung is consistent with the color dipole 
picture in the limit where the number of colors N c is large. The color dipoles are quark-antiquark 
pairs which can radiate arbitrarily many soft gluons, and the evolution consists in the splitting 
of any such a dipole into two. We construct the color glass weight function of an onium as a 
superposition of color dipoles, each represented by a pair of Wilson lines. We show that the action 
of the Bremsstrahlung Hamiltonian on this weight function and in the large- A^ c limit generates 
the evolution expected from the dipole picture. We construct the dipole number operator in the 
Hamiltonian theory and deduce the evolution equations for the dipole densities, which are again 
consistent with the dipole picture. We argue that the Bremsstrahlung effects beyond two gluon 
emission per dipole are irrelevant for the calculation of scattering amplitudes at high energy. 
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1 Introduction 



Our purpose in this paper is to demonstrate that the perturbative evolution of a dilute 
hadronic system with increasing energy can be effectively described as the evolution of a 
system of color dipoles provided we consider the limit in which the number of colors N c 
is large and that the energy remains low enough for the evolved system to be still dilute. 
By 'dilute' we mean that the hadronic system is non-saturated : the gluon density is low 
enough for the recombination processes to be unimportant. 

At an abstract, wavefunction, level, where the evolution is viewed in terms of the partonic 
content of the system and not of its interactions, the effectiveness of the dipole picture 
has been of course demonstrated in the original paper by Mueller [1]. Here, however, we 
are primarily interested in the problem of scattering — the dilute hadronic system is the 
target which scatters with an external projectile, itself assumed to be dilute — and to that 
purpose we need to specify not only the evolution of the partons in the target, but also 
the way how an individual parton (gluon or dipole) couples to the projectile. In previous 
applications of the dipole picture to scattering [2-7] , one has always assumed that a target 
dipole which partakes in the collision exchanges exactly two gluons — i.e., it undergoes 
single scattering — with the projectile. Whereas this approximation is indeed justified 
at relatively high energy, where the large dipole density in the target favors the multiple 
scattering with different dipoles [2] (see also the discussion in Sect. 5 below), this is not 
really correct at low energies, where e.g. the double scattering off a same target dipole via 
four gluon exchange competes with the scattering off two different dipoles, via twice two 
gluon exchange. 

Notice that multiple gluon exchanges between a target dipole and the projectile correspond 
to many gluon radiation from the quark and the antiquark legs of the dipole, that is, to 
gluon Bremsstrahlung. The theoretical description of the Bremsstrahlung in the high- 
energy evolution of a dilute QCD system has recently became available [8-10], but its 
relation with the dipole picture at large N c has not been established so far, in spite of 
previous attempts [11] which emphasized the complexity of the problem. It is our present 
objective to fully clarify this problem by showing that, at large N c , the evolution generated 
by the Bremsstrahlung Hamiltonian in Refs. [8-10] can be recast in the language of the 
dipole picture. 

Specifically, we shall find that, for large N c , the wavefunction of the dilute target (the 
'onium') can be described as a collection of color dipoles, where each dipole is allowed to 
radiate arbitrarily many small-x gluons in the eikonal approximation. At a mathematical 
level, a dipole is represented by a pair of Wilson lines (one for the quark and the other 
one for the antiquark) built with the color field of the radiated gluons. The evolution 
of the system with increasing energy is governed by the action of the Bremsstrahlung 
Hamiltonian on the 'onium' wavefunction. As we shall see, at large N c this evolution 
proceeds through dipole splitting: one pair of Wilson lines splits into two such pairs which 
have one common transverse coordinate (at the position of the emitted gluon). 
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The formulation of the dipole picture that we shall naturally arrive at is of the color 
glass type [12-14] : The dilute target is effectively represented as a stochastic ensemble of 
classical color fields — the fields radiated by the dipoles — which are distributed according 
to a functional 'weight function' whose evolution we shall compute. This evolution can 
be reformulated in terms of probabilities for the dipole configurations (or, equivalently, 
in terms of dipole n-body densities, with n > 1), and then it reduces to the original 
wavefunction evolution 2 by Mueller [1], as expected. But the color glass formulation 
turns out to be more convenient for applications to the scattering problem, as it gives 
directly the distribution of the target color fields to which couple the projectile. 

At this point we should remind that a color glass formulation of the dipole picture has 
been already given in Refs. [3,5], but under the assumption that each dipole can radiate 
only two gluons. When performing the corresponding approximation on our subsequent 
results (this amounts to expanding the Wilson lines to lowest non-trivial order), we shall 
recover the weight function of Refs. [3,5], as expected. However, this does not mean that 
the Bremsstrahlung Hamiltonian -£/brem should mechanically reduce at large N c to the 
'Dipole Model' Hamiltonian i^Msw introduced by Mueller, Shoshi and Wong [5]. Rather, 
these two Hamiltonians describe the same physical process — the splitting of one dipole 
into two — but by acting on different dipole operators: -£/msw acts m the Hilbert space 
of bare dipoles, by which we mean the dipoles which are allowed to radiate only two 
gluons, whereas -^brem acts on dipoles which are fully dressed by the radiation. In fact, 
the "large-iV c limit of -£/brem" is not a well-defined concept by itself: The simplifications 
appropriate at large N c can be performed only in the process of acting with H^kem on 
the onium weight function, or on gauge-invariant correlations expressing dipole densities. 

Since both Herem and Hyisw describe dipole splitting (at large N c ), they are equivalent in 
so far as the evolution of the onium wavefunction is concerned: they generate the same evo- 
lution equations for the dipole densities (or the associated probabilities), which moreover 
coincide with the corresponding equations previously derived within the 'abstract' dipole 
picture [1-4, 16]. However, since the physical information encoded in the dipole operator 
is different in the two theories — this is richer for -£^brem, where a single dipole is allowed 
to radiate arbitrarily many gluons — , and moreover this information is relevant for the 
scattering, it follows that the evolution equations for scattering amplitudes generated by 
the two Hamiltonians will be different in general. For instance, the scattering between the 
onium and two external dipoles may proceed either via twice two gluon exchange with 
two different internal dipoles, or via four gluon exchange with a single internal dipole. 
Hmsw describes only the former process, and thus generates the 'fluctuation' terms in 
the evolution equations with Pomeron loops [4,6], while -^brem describes both of them, 
and thus in principle it generates more complete equations. Still, as mentioned before, 
and will be demonstrated in Sect. 5 below, the latter process is suppressed at high energy 
with respect to the former one, and thus is irrelevant for the study of the approach to- 
wards saturation and the unitarity limit. A similar conclusion has been recently reached 



2 More precisely, the formulation of the 'onium' wavefunction evolution which will naturally 
emerge from our calculations is that in Refs. [3,4, 15, 16]. 
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by Marquet, Mueller, Shoshi and Wong [17, 18]. 



The relation between -^brem and Hmsw alluded to above is in fact similar to that between 
the JIMWLK Hamiltonian and its 'Pomeron' approximation at large N c , as discussed in 
Ref. [7]. The JIMWLK Hamiltonian [13,19,20] describes non-linear effects (gluon re- 
combination, or merging) in the evolution of the target wavefunction towards saturation. 
Alternatively, when acting on scattering operators built with Wilson lines, it describes 
gluon splitting in the projectile followed by the scattering between the products of this 
splitting and the target. At large N c , it is convenient to consider a projectile which is 
itself made with dipoles. Then, as shown in Ref. [7], the action of -£/jimwlk on a dressed 
external dipole (where the 'dressing' now refers to multiple scattering, as encoded in the 
Wilson lines) is equivalent — in the sense of generating the same evolution equations for 
the scattering amplitudes at large N c — to that of a simpler 'Pomeron merging' Hamil- 
tonian on a bare external dipole, i.e., a dipole which exchanges only two gluons with the 
target. This close correspondence between the two problems — the JIMWLK evolution in 
the high density regime and the evolution including Bremsstrahlung in the dilute regime 
- should not came as a surprise, in view of the duality between .Hbrem and -£/jimwlk on 
one side [9, 10], and between i^Msw and the Pomeron-merging Hamiltonian on the other 
side [7]. Our subsequent analysis will shed more light on this correspondence, by show- 
ing that the duality holds, more precisely, between the action of Hjimwlk on scattering 
operators for external dipoles and, respectively, that of -/^brem on radiation operators for 
the internal dipoles. In brief, the action of -^jimwlk on the projectile is dual to that of 
-^brem on the target. 

The plan of the paper is as follows: In Sect. 2 we briefly review the Hamiltonian formulation 
of the high-energy evolution in the dilute regime and in the presence of Bremsstrahlung, 
then explain the action of -?/brem in terms of Poisson brackets (or 'commutators'), and 
finally demonstrate as a simple exercise that the action of this Hamiltonian on the 2-point 
function of the color charge density generates the BFKL equation, as expected. Sect. 3 
contains our main results, namely the construction of the onium weight function in terms 
of dressed dipoles, and the proof of the fact that the action of .Hbrem on this weight 
function generates the dipole evolution at large N c . In Sect. 4, we show how to express 
dipole number densities in terms of gauge-invariant correlations of the color charge. We 
introduce the dipole number operator, and verify that, under the action of .Hbrem (at large 
N c ), the dipole densities obey the evolution equations expected in the dipole picture. In 
Sect. 5 we consider the onium scattering with two external dipoles and show that the 
elementary processes which involve the double scattering of a same target dipole are 
subdominant at high energy. Finally, Sect. 6 contains our conclusions. 



2 QCD evolution in the low density regime: Bremsstrahlung Hamiltonian 

The physical problem that we have in mind is that of the scattering between a dilute 
hadronic system propagating in the positive z (or positive x + ) direction — the right- 



4 



moving target — and a system of color dipoles moving in the negative z (or positive x~) 
direction — the left-moving projectile. When the projectile is made of a single dipole, the 
scattering measures the gluon distribution in the target, i.e., the average gluon number 
density. A projectile made with several dipoles will also probe fluctuations in the gluon 
number. In the limit where the number of colors N c is large, we shall be able to effectively 
describe the gluon distribution in the target in terms of (internal) dipoles. Then the 
scattering with the external system of dipoles will measure the dipole distribution in the 
target, that is, the average dipole number density and the corresponding fluctuations. 

The operator expressing the scattering amplitude for an external dipole (x, y) in the 
eikonal approximation reads: 



T(x, y) = l-± r tr(v\x)V(y)) ~ ^ (a a (x) - a a {y))\ (2.1) 



where the second, approximate, equality holds in the case where the target is dilute, which 
is the case of interest for us here. In this equation, the Wilson line 

V\x) = P exp ji0 J dx~a a (x-, x) t a | , (2.2) 

(the t a 's are the generators of the SU(iV c ) algebra in the fundamental representation and 
the symbol P denotes path-ordering in x~) describes the eikonal scattering between a 
quark with transverse coordinate x and the light-cone component A\ ) of the 

color field in the target. Similarly, V(y) describes the scattering of the antiquark in the 
dipole. Furthermore, 

a a (x) = J dx~a a (x~,x) (2.3) 

is the effective color field in the transverse plane, as obtained after integrating over the 
longitudinal profile of the target, and is related to the corresponding color charge density 
in the target p a (x) via the two-dimensional Poisson equation 

-Via o (a0 = Pa (x) . (2.4) 

Note that the time variable x + is suppressed in the equations above, since the scattering is 
quasi-instantaneous; thus, a and p are the color field and charge at the time of scattering, 
which is x + = oo in the conventions of Ref. [10] : a = and p = p^. 

The physical scattering amplitude for a projectile made with k dipoles is obtained as 
(T^) T = (T(xi,yi)T(x2,y2) • • -T(a;i,j/fc)) r , where r ~ Ins is the rapidity gap between 
the projectile and the target (note that we use a Lorentz frame in which most of the 
total energy is carried by the target) and the brackets denote the average over the target 
wavefunction, which in the spirit of the color glass formalism [14] is computed as an 
average over the color charge p with weight function Z T [p\. E.g., for a single dipole: 

(T(a!,y)> T ~ / D[p] ^ (a a (x) - a a (y)) 2 Z T [p] ■ (2-5) 
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To account for the correlations induced by Bremsstrahlung, the color charge density should 
be treated as a time-dependent variable [8-10], so Z T [p] is a functional of p a (x + , x) (see 
below for more details). 

According to the above equations, in order to compute physical scattering amplitudes and 
their evolution with r, it is sufficient to know the target weight function Z T [p] and the 
corresponding evolution equation. The latter can be compactly written as 

d 

-Q^ z r[p] = - H BREM Z T [p], (2.6) 

with the following Hamiltonian for color glass evolution in the dilute regime and in the 
presence of Bremsstrahlung [8-10] : 



xyz 



(For more clarity, we have temporarily restored the subscript oo denoting the x + variable 
of the field p.) In this equation, 



K(x,y,z) 



z ) ■ [y 



x 



y) 



2 ' 



(2.8) 



Furthermore, W and are all-order differential operators defined as, e.g., 

W(x) — T exp \—g I dx+T a 6 ) (2.9) 
{ J 5p a {x+,x)) 

where T denotes time-ordering and the color matrices T a are in the adjoint representation. 
Physically, W and W' are Wilson lines describing the radiation of an arbitrary number 
of small-x gluons in the eikonal approximation from a single gluon with a larger value 
of x. (As usual, x denotes the longitudinal momentum fraction carried by a gluon, with 
r = ml/a;; see Ref. [10] for details.) When Eq. (2.9) is expanded in powers of 5/5p, each 
such a power describes the creation of a source for the emission of one gluon (cf. Eq. (2.4)). 
Since these sources are dilute, the radiated fields are never strong. However, by keeping 
terms of all orders in S/Sp in the Wilson lines, and thus in the Hamiltonian (2.7), one 
includes in the evolution processes like 2 — > n gluon splitting through which the n-point 
functions of p with n > 2 get built from the 2-point function (pp) in the dilute regime. 

For what follows, the explicit representation (2.9) for W is not so important. Rather, what 
matters is the structure of the Poisson brackets which define the action of the Hamiltonian 
(2.7) on the Hilbert space of the operators 3 Ofpoo, W\. 

3 In what follows, we shall use the notation W for the temporal Wilson line in a generic, or in 
the fundamental, representations, and we shall keep the more specific notation W for the adjoint 
representation. 
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p^x), W bc (y)} = g(T a W(x)) bc 5^(x - y), 
p^x), pKy)] = -igfto&ix) 5^\x - y), 
W ab (x),W cd (y)} = 0. (2.10) 

The first two equations above show that the color charges act as infinitesimal gauge 
rotations of the Wilson line W at its end points, that is, as Lie derivatives. One can 
check that these commutation relations satisfy the right properties expected for Poisson 
brackets, in particular, they obey the Jacobi identity. 

The non-commutativity of the color charge variables p^x) with themselves is a source of 
potential difficulties (like ambiguities in the ordering of the operators) in the construction 
of the correlation functions and of the color glass weight function [8-11]. Still, as we shall 
see, this problem can be systematically avoided in the large-iV c limit, where we shall be 
able to explicitly construct the weight function Z T [p] in a way which is free of ambiguities. 

The formal evolution equation for some operator Ofp^l/F] is defined by its Poisson 
bracket with H BREM : 

^-0[ PoD ,W] = [H BREM ,0]. (2.11) 

Note that for a given ordering of the operators which compose O, the operation above is 
unambiguous. However, its result will generally change when we permute the operators 
within the definition of O. Once again, this difficulty will not show up at large N c . 

One can check that, for physical observables at least, Eq. (2.11) is indeed consistent with 
the equation (2.6) for the evolution of the weight function together with the definition 
(2.5) of the color glass average. To see this, note first that the physical observables are 
gauge-invariant, which in the present context means that they are invariant under the 
gauge transformations dependent upon x + [10]. When the action of -£Tbrem is restricted 
to such gauge-invariant operators, Eq. (2.7) can be equivalently replaced by [21] 

#brem = / M xyz [l + W x W y - W X W Z - W z W y ] ab PloWpKy), (2.12) 

xyz 

with M xyz denoting the dipole kernel [1] : 

(x - y) 2 

M(x, y, z) = _ = K, xxz + JC yyz - 2IC xyz . (2.13) 

[x z) [z y) 

(This can be proven through the 'dual' version of the arguments used in Ref. [21] for the 
case of the JIMWLK Hamiltonian.) From Eq. (2.10) one can check that, in the presence of 
the dipole kernel, the color charge operators in Eq. (2.12) can be freely commuted through 
the Wilson lines there; in writing Eq. (2.12) we have used this freedom to commute both 
factors of p fully to the right, which is convenient for the subsequent manipulations. Also, 
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the non-commutativity between p^o(cc) and p b 00 {y) plays no role in Eq. (2.12) since the 
Wilson-line part of the integrand, 



ah 



h ab (x, y, z)=[l + W x Wl - W x Wl - W z Wl 
is symmetric under the simultaneous exchange a <-> b and x <-> y. 



(2.14) 



Note furthermore that, if the operator OI/^W] contains a factor of W on the left of 
all the other operators, then when computing the color glass average (0[poo, W]) T this 
factor can be replaced by one. Indeed, when expanding the exponential in Eq. (2.9), all 
the terms but the first one yield total derivatives which vanish after integration over p. In 
the subsequent manipulations, it will be often convenient to 'normal-order' the operators 
by pushing factors of W all the way to the left and the replacing them by one when 
computing the average. 

We are now prepared to check that Eqs. (2.6) and (2.11) are consistent with each other. 
To that aim, take the color glass average in Eq. (2.11). This involves 



h ab (x,y,z)p a 00 (x)p b 00 (y), 0\ Z T [p\ 
J D[p] (K^plO - Oh^pl) Z t [ P ] 
J D[p] ( - Oh^pl) Z T [ P ] — > / D[p] O (- #brem^[pD , 



(2.15) 



where we have used the fact that h[W] — > when W — > 1. By the same argument, we 
deduce that the average of Eq. (2.11) involves only the commutator [h, O] : 



-1 
16^3 



D[p] / M xyz h ab (x,y,z),0 p^p^y) Z T \p\. 



(2.16) 



xyz 



As a first, relatively simple, application of the above formalism, let us derive in this way 
the BFKL equation. We shall argue later that, at large N c and for x ^ y, the charge- 
charge correlator (p a (x) p a (y)) T is proportional to the dipole number density n T (x,y), to 
be precisely defined in Sect. 4. Here and in what follows, we use the simpler notation 
Px = Pto( x )- Specifically : 



{p a {x)p a (y)) T = -g 2 C F n T (x,y) +n T (y,x) for x^y. 



(2.17) 



Since, on the other hand, the BFKL equation for a 2-point function is well known to 
emerge independently of the large- N c approximation, we expect the quantity in Eq. (2.17) 
to obey the BFKL equation for arbitrary N c . Let us check that this is indeed the case. 
According to Eq. (2.16), we have: 



M 

uvz 

( \ h cd {u,v,z), 

PxPy 



Pupt 



(2.18) 
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where we have relabeled the transverse coordinates internal to .Hbrem as u, v and z. By 
repeated use of the commutation relation, we shall move h c d (which contains the Wilson 
lines W) to the left and then set W = 1 



[h cd ,p a x p a y ] = [h cd ,p a x ]p a y + p a x [h cd ,p a y ] 

= ^Cd, plVy + [hcd, Py]p a x + [Pi, [hcd, Py] ] 
— [h Cd ,Py]], 



(2.19) 



where, as shown in the last line, only the double commutator must be retained when 
computing the color glass average in Eq. (2.18). Indeed 



[h cd ,p a x \ 



= 0, 



(2.20) 



as it can be easily checked by using Eq. (2.14) together with the following commutator, 
which in turn follows from Eq. (2.10) (S ux = S^(u — x)) : 

[p%, [WjvlU] = g5 ux (T a W u Wi) cd - g5 vx (W u WlT a ) cd . (2.21) 

To evaluate the right hand side of Eq. (2.18) we also need 

Pi, [P y , {W u Wl) cd ]\ =g 2 5 uy {5 ux (T b T a W u Wl) cd - 5 vx (T b W u WlT a ) J 

-g 2 5 vy {5 ux {T a W u W v T b ) cd ~5 vx {w u WlT a T b )J . (2.22) 

From now on, simple algebra yields 



d a 2 N 



J M um ( - W„„ + KJ,y + KAy) (f&Vlr , (2.23) 



where we have neglected the terms which vanish at x ^ y and used the symmetry u <-> v 
of the kernel in the Hamiltonian. After also using Eq. (2.17), this is finally rewritten as 



dn T (x,y) _ 
Or ~ 2tt 



-M(x,y,z) n T (x,y) 

+ M(x, z, y) n T (x, z) + M(z, y, x) n T (z, y) 



(2.24) 



(with a s = a s N c /ir), which at large N c is recognized as the BFKL equation for the dipole 
number density [4, 16], but which is valid as written for arbitrary N c . 

Note that the same equation could have been obtained by first expanding the Wilson lines 
in the Hamiltonian (2.7) in a power series in derivatives, then keeping the first non-trivial 
terms (the second order ones) in this expansion to deduce the BFKL Hamiltonian: 



BFKL — 



167T 3 



/ 



M uvz f ace f b 



5p c (u) Sp c (z) 



5p d (z) 5p d (v) 



p a (u)p b (v), 
(2.25) 
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and finally using this Hamiltonian in the evolution equation for (p a (x)p a (y)) T : 

— (p«( x )p»(y)) T = J D[p]p a ( X )p a (y)(- H BFKL Z T [p]) • (2-26) 

After some integration by parts, the functional derivatives in Eq. (2.25) are brought to 
act on the factors of p in the operator, and then Eq. (2.23) immediately follows. 

One may think that the first derivation of Eq. (2.23), in which the Wilson lines inside 
-^brem were kept unexpanded, is more general than the one based on the BFKL Hamilto- 
nian (2.25), but this is only illusory: Since O = p a x p a y is quadratic in p, all the higher order 
p-derivatives beyond the second-order ones kept in Eq. (2.25) do not contribute to its 
evolution. Thus, in so far as Eq. (2.23) is concerned, the two methods presented above - 
the use of the commutation relations and the derivative expansion of the BREM Hamil- 
tonian — are equivalent with each other. However, for more general situations (e.g., more 
complicated observables, or the evolution of the weight function that we shall consider in 
the next section), the use of the commutation relations turns out to be more convenient 
as it avoids potential ambiguities with the x + -ordering of the functional derivatives in 
the expansion of the Wilson line (2.9). 



3 Dipole picture from the Bremsstrahlung Hamiltonian 

In this section we shall derive the dipole picture from the action of the Bremsstrahlung 
Hamiltonian in the large-iV c limit. More precisely, we shall show that, if one starts with 
a single color dipole — a quark-ant iquark pair with the quark at u and the antiquark 
at Vq — at the initial rapidity r = 0, then the partonic system produced at some higher 
rapidity r through the evolution described by H BREM at large N c can be itself characterized 
as a collection of qq color dipoles (an 'onium'), which evolves through dipole splitting. 
Besides the Bremsstrahlung Hamiltonian (2.7), the crucial ingredient in this picture is the 
description of the onium wavefunction as a color glass, which in turn requires the proper 
definition of the dipole operator as a color source. 

The first formulation of the onium as a color glass has been given in Ref . [3] , under the 
assumption that each dipole is a color source for only two gluons. The evolution of the 
associated weight function through gluon splitting has been exhibited in Ref. [3], but it was 
only later, in Ref. [4], that one has realized that this evolution cannot be fully accounted 
for by the JIMWLK Hamiltonian. The appropriate Hamiltonian has been constructed 
shortly after, by Mueller, Shoshi and Wong [5]. It is the sum of the BFKL Hamiltonian 
(2.25) plus a term involving four p-derivatives which describes dipole splitting (again, 
under the assumption that each dipole can radiate only two gluons). In what follows, 
we shall generalize the construction in Refs. [3, 5] to the case where a dipole can radiate 
arbitrarily many (small-rr) gluons. That is, we shall construct the corresponding dipole 
operator and onium weight function, and show that the Bremsstrahlung Hamiltonian (2.7) 
is the appropriate generalization of the MSW Hamiltonian [5]. 
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We start with a brief summary of the results in Refs. [3, 5], to which we shall refer as the 
'Dipole Model' (DM). The onium (color-glass) weight function in the DM reads: 



oo „ N 

Z? M [P] = E dT N P N ({z i };r)l[Dl(z i _ 1 ,z i )5[p], (3.1) 

N=l J i=l 



where P]y({zi};T) denotes the probability density to find a given configuration of N 
dipoles at rapidity r (the configuration being specified by N — 1 transverse coordi- 
nates {zi} = {zi, z 2 , ...Z/v-i}, such that the coordinates of the N dipoles are (z ,Zi), 
(zi, z 2 ),...,(z N _i, z N ), with z = u and z N = v Q ), and dT N denotes the measure for the 
phase-space integration: dT N = d 2 z 1 d 2 z 2 . . . d 2 z N _i. Furthermore, D\{u,v) is the DM 
dipole creation operator, with the dipole assimilated to the source of two gluons: 

i&», y ) S i + £( * - * y, 0.2) 



AN c \5p a (x) 5p a (y) 

where p a (x) = J dx~p a (x~,x) should be interpreted as the color charge density in the 
transverse plane at the interaction time. (There is no explicit x + -dependence in the DM 
picture; see also Eq. (3.9) below.) Finally, the delta-functional S[p] = 5[p a (x)] is defined 
in the context of the functional integral over p a (x). 

The evolution of the DM is driven by 

^Z™[p] = -H MSW Z™[p], (3.3) 



with the Mueller-Shoshi-Wong Hamiltonian 



Hmsw = J M{x,y,z) -Dl(x,y) + Dl(x,z)Dl(z,y) 



D (x,y), (3.4) 



xyz 



where D (x,y) is the dipole annihilation operator within the DM: 



D (x,y) = -— p a (x)p a (y) for x + y. (3.5) 



1 

One can check indeed that: 

[D (x,y),Dl(u,v)} w i (S ux 8 vy + 5 uy 8 vx ) , (3.6) 

where the approximate equality sign means that the equality holds in the large-iV c limit. 
(In Sect. 4, we shall demonstrate a relation similar to Eq. (3.6) in a more general context.) 

By inserting Eqs. (3.1) and (3.4) into Eq. (3.3), then using Eq. (3.6) to successively 
commute the annihilation operator D (x,y) from Husw to the right of the creation 
operators Z?o(zj_i, zi) from Z® M , and finally using the fact that Dq5[p] = 0, one finds the 
following evolution equation for the onium weight function: 
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c\ — oo N 

dr 271 N=l J i=l { 



x [£>S(«i-i, *) - ife-i, *)] II ^fe-i, *j) S\p], (3.7) 

which is the expression of the dipole picture in the color glass representation [3]: the 
evolution of the weight function proceeds through dipole splitting. From here on, one 
can proceed in the standard way [3, 4, 16] to deduce evolution equations for the dipole 
densities and probabilities (see the discussion at the end of this section and in Sect. 4). 

We now return to the general case where a dipole can emit an arbitrary number of gluons, 
and show that an evolution equation similar to Eq. (3.7) (which is synonymous of the 
dipole picture) is obtained also in that case provided we replace Husw by .Hbrem and the 
'bare dipole' creation operator D\ by the following, 'dressed dipole', operator [11] : 

D\x,y) = j r tr(w(x)W\y)), (3.8) 

where the Wilson lines are in the fundamental representation (W stays for the quark, and 
for the antiquark). To second order in the expansion of the Wilson lines in powers of 
S/Sp, and with the following identification 



5p a (x) J 5p a (x + ,x) 

Eq. (3.8) reduces to the two-gluon emission operator, Eq. (3.2). In general, Eq. (3.8) 
can be seen as the gauge-invariant generalization of Eq. (3.2) to the regime where the 
derivatives are formally strong : g J dx + (5/5p) ~ 1 (cf. the discussion after Eq. (2.9)). 

The onium weight function will be now constructed by analogy with Eq. (3.1). That is, 
we start by assuming that the weight function Z T [p] in the large- N c limit can be cast in 
the following, dipolar, form : 

oo „ N 

zap] = E dr N p N ({z i y,T)]lD\z i . 1 ,z i )s\p], (3.10) 

JV=1 i=l 

and then show that this particular structure is indeed preserved by the evolution under 
the action of the Bremsstrahlung Hamiltonian (2.7) and for large iV c . 

The subsequent mathematical manipulations can perhaps be better understood if one no- 
tices their analogy (in fact, duality) to manipulations which are by now familiar in the con- 
text of the JIMWLK evolution (see, e.g., Refs. [13,21]). Specifically, the Bremsstrahlung 
and JIMWLK evolutions are known to be dual to each other [9, 10], in the sense that the 
corresponding Hamiltonians and also the respective Poisson brackets get interchanged 
with each other under the following duality transformations 4 : 



Note that the "oo" subscript in refers to x , and not to x + ; see Ref. [10] for details. 
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hafej^AW. Vt(x) ~ W(X) 



(3.11) 



Similarly, the dipole creation operator (3.8) is dual to the S'-matrix operator S(x, y) = 
tr(V^(x)V(y)) which describes the scattering of an external dipole in the high-density 
regime specific to the JIMWLK evolution (cf. Eq. (2.1)). Thus, clearly, -^brem acts on the 
internal dipoles in the same way as -£/jimwlk does on the external ones. We know already 
that the action of -^jimwlk on scattering operators generates the Balitsky equations [22] . 
At large N c and for projectiles built with dipoles, the Balitsky equations close in the 
space of dipole operators, and are consistent with the dipole picture for the projectile 
wavefunction [4, 16,23]. This lets us anticipate that the action of -£^brem on the dipolar 
weight function (3.10) at large N c should similarly generate the dipole picture for the 
target wavefunction. This will be verified explicitly in what follows. 

We thus need to evaluate 



t^T z Ap] = - H BREM Z T [p] 

= £ J dTNPN ^Y^-3 J M^K.pIpIUD^z^, zt) 5[p] (3.12) 

N=l 71 xyz i=l 

in the large-iV c limit. For more clarity, it is convenient to use the notations Ui = z%-\ 
and Vi = Zi for the transverse coordinates of the i-th dipole; it is then understood that 
u i+1 = Vi. Let us focus on the action of H BREM on the weight function Z N for a given 
configuration of N dipoles. This involves: 



N 



h abP a x p y Z N ({ Ui ,V t }) EE hatplPyHD^U^vJSip}. 



(3.13) 



i=l 



To proceed, we need to commute p x p y to the right of the dipole creation operators and 
then use pS[p\ = 0. The relevant commutator is 



o a o b D ] 

HxHy) 



(f £)t ] J> , [ b D \ 1 a , [ a l b £>f 
rxi uv I y 1 ry ^ uv rx 1 I xi ryi ^ v. 



By adapting Eqs. (2.21)-(2.22) to the fundamental representation, we obtain 



o a o b L> f 

rxryi ^ uv 



-- $ X vMt a W u Wl) p\ 

+j r (S yu - 8 yv )tr{t b W u Wl) pi 



+g 2 ($y U - Syv) 



S xu ^tv(W u Wlt b t a ) - 5 xv ±-tr(W u Wlt a t b ) 



(3.14) 



(3.15) 



where the expression in the last line represents the double commutator term. We shall 
shortly verify that the first two terms on the right hand side (proportional to p), when 
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multiplied by h a b and acting on Z N , give subleading contributions at large N c compared 
to the last term, which does not involve p. Introducing the more compact notation 



o a o b -D f 



. , = i^(5 xu -5 xv Mt a W u Wl) J > 

non— dipole _/V c 

9 



Py 



+f- c (Sy U -S yv )tr(t b W u Wl)p a x 
= A a rfi 4- R b n a 

uvx ry ' uvy rxi 



(3.16) 



where A and B do not contain p, we can write the result of the first commutation as 
follows (in simplified notations whose meaning should be obvious) : 



n a n b 71* 71* D* — n a n b D* 71* -I- Tn a n 6 71* 
PxPy U \ U 2 ■ ■ ■ U N — U \ PxPy U 2 • • • U N + PicPj/i ^Atiwi 



non— dipole 



H* .H* 



N 



+ 



Pi, 



p b L>* 

ryi u\ vi 



£>!...£>t. 



(3.17) 



The last, double-commutator, term in the r.h.s. of the above equation is already indepen- 
dent of p, so we can operate with h a b on it : 



1 



167T 3 



M xyz h ab (x,y,z) p a x , p b y ,Dl 



UlVl 



xyz 



■g^- / M UlVlZ h ab (u 1 ,v 1 ,z)tY(W Ul Wl l t b t a ) 



= ^ / M U1V1Z { - ^(W^Wi) + ^(W Ul W z ) ^tr^W*)} (3.18) 



where we have also used the identity {t a )ij(t a ) k i = \8u8jk — -^fS^Ski- The last equation is 
dual to the first Balitsky equation [22], and in fact this has been obtained here through 
manipulations similar to those usually performed in the derivation of the Balitsky hier- 
archy from the JIMWLK equation [13]. The first, negative, term within the braces in 
Eq. (3.18) describes the probability that the original dipole (u±, Vi) survive without split- 
ting, while the second, positive, term describes the splitting of the original dipole into the 
new dipoles {ui,z) and (z,Vi) . 

We shall now check that the non-dipolar contribution to Eq. (3.17) is indeed suppressed 
at large N c . To that aim, we take one particular piece in the non-dipolar commutator 
(3.16), say the first piece A^ iViX Py, and consider its action on the dipole creation operators 
which appear on its right in Eq. (3.17): 



A a n b D* /I* /I* r)\n\ = /I* A a n b /I* /I* r)\n\ 

uivix ry U2V2 U3V3 " " " u^vm irl U2V2 u\ v\ x ry U3V3 • • • umvn Lri 



+ 



At 



p b L>* L>* 

u,\V\x I yi U2V2 



U3V3 



Di NV J[p}. (3.19) 



The above commutator is independent of p, and the same is true for all the other commu- 
tators generated when p b is further commuted towards the right. It is therefore sufficient 
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to evaluate the action of h ab on one such a commutator. We have: 



A a UlVlxP b y ,Dl 2V2 = -^(5 XUl - 5 XVl )Tr{t a W Ul Wl) P h y ,D{ 

2 



!J (<W - S xvl )(5 yU2 - 5 yV2 )Tr(t a W Ul Wl)Ti(t b W U2 Wl)- (3-20) 



c 



Let us consider the action of the first two terms in h ab (x, y, z) on this commutator: 



xyz 



M xyz (5 XUl - 5 XVl )(5 yU2 - 5 yV2 )(l + W x Wl) a Ti(t a W Ul WljTi(t b W U2 Wi 2 ) 



2 

= j- I (M UlU2Z - M UlV2Z - M VlU2Z + M V1V2Z ) ^Ttr(W Ul WlW Ul Wl 2 ) 



2N C 

+^(Wl l W ul Wl 2 W U2 ) - ^tr(W M1 Wi) ±-tr(W U2 Wl 2 : 

(3.21) 

The r.h.s. of Eq. (3.21) is of order a s /N^, and hence it is suppressed by a factor of 
1/iVg compared to the dipole contribution in Eq. (3.18). Similarly, one can verify that 
the other contributions, due to the last two terms in h a b(x,y, z), are also suppressed. 
Once again, this property has a dual counterpart in the context of the Balitsky-JIMWLK 
equations: The non-dipolar terms in the evolution equation for the scattering amplitude 
of a projectile made with two dipoles are suppressed at large N c . 

Retaining only the dipolar contribution in Eq. (3.18), we finally obtain the following 
evolution equation for the weight function in the large-iV c limit: 



<9 a 00 r N r 

Z T [p}=7rE dT N P N (r)J2 M(u, t ,v t ,z) 

dr 2n N=l J i=l { 



x [D\u u z) D\z, Vi ) - D\u h Vi )\ n D\ Uj , Vj ) 5[p}. (3.22) 

As anticipated, this has the same structure as the corresponding equation in the Dipole 
Model, Eq. (3.7), except for the replacement of the creation operator for a bare dipole D\ 
with the corresponding operator for a dressed dipole . Note that the non-commutativity 
of the color charge operators p a (cf. the second equation (2.10)) did not play any role in 
the manipulations leading to Eq. (3.22). This suggests that the ordering of the operators 
should be irrelevant in the large- N c limit (at least, for the dipole-related variables). This 
conclusion will be further supported by the developments in Sect. 4. 

From Eq. (3.22), the dipole picture of the target wavefunction can be developed along the 
same lines as in Ref. [3]. Namely, Eq. (3.22) is consistent with the dipolar structure of the 
weight function, Eq. (3.10), provided the probability densities Pn which enter the latter 
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obey the following Master equation : 



dP/vQi, ...zjy-i; t) _ Os_ 
Or ~ 2tt 



N 

/ d 2 z M{zi^Zi,z) 
i=i J 



P N (z 1 , ...z N -i,t) 



a, 



+ tt~ Y] M(zi-i, z i+ i, Zi) P N _ 1 (z 1 , Zi_i, z i+1 , z N -i,t). 
2n i= i 

(3.23) 

The first term in the r.h.s., proportional to P/v, is the loss term, which describes the 
splitting of one dipole from the original configuration of N dipoles. The other terms, 
proportional to Pjv-i, are Q^in terms showing the formation of the iV-dipole configuration 
of interest through the splitting of one dipole in original configurations of N — 1 dipoles. 

One can summarize the previous discussion as follows: At large N c , the Bremsstrahlung 
Hamiltonian (2.7) acts on the onium weight function built with dressed dipoles, Eq. (3.10), 
in the same way as the Dipole Model Hamiltonian Pmsw ac ts on the weight function (3.1) 
built with bare dipoles. But this does not imply that Pmsw can be deduced from Pbrem 
by somehow taking the large-iV c limit of the latter. That is, the derivative expansion 
of Hbrem, Eq. (2.12), up to fourth order in S/Sp is neither equivalent to the expression 
(3.4) of Hmsw, nor it reduces to it in some suitable large-iV c limit. In fact, since the 
expressions for the onium weight function in the Dipole Model, Eq. (3.1), and in the 
general Bremsstrahlung case, Eq. (3.10), differ from each other by terms of order (5/5p) 3 
or higher, the respective Hamiltonians must differ from each other too (already at order 
(5/5p) 3 ), in order to yield identical results when acting on different dipole operators. This 
observation explains, in particular, why the analysis in Ref. [11], which has limited itself 
to the extraction of the four derivative terms in the Bremsstrahlung Hamiltonian (2.7), 
has met with difficulties when trying to generate the dipole picture from Pbrem- 



The above discussion finds its dual counterpart in the analysis of the JIMWLK evolution 
in Ref. [7] : At large N c , Pjimwlk acts on the dipole operator S(x, y) — jj- tr(V^(x)V(y)), 
which describes multiple scattering, in the same way as the 'Pomeron' Hamiltonian ob- 
tained [7] as the dual partner of the MSW Hamiltonian (3.4) acts on the 'bare' operator 
So(x,y) = 1 — (a a {x) — a a (y)) 2 , which describes single scattering via two gluon ex- 
change (cf. Eq. (2.1)). The 'Pomeron' Hamiltonian at high density and, respectively, the 
Dipole Model Hamiltonian (3.4) in the dilute regime are just effective Hamiltonians, which 
generate the correct evolution equations at large N c (for dipole scattering amplitudes in 
the first case, and for target dipole densities in the latter), but by acting in a simplified 
Hilbert space in which each dipole is characterized by the exchange of only two gluons. 

It is finally interesting to notice the different ways how the actual QCD Hamiltoni- 
ans (BREM and JIMWLK) and their effective counterparts (MSW and, respectively, 
Pomeron) act in the corresponding Hilbert spaces. By construction, Husw implements 
the dipole evolution in the most straightforward way: when acting on a collection of (bare) 
dipoles, it first annihilates one dipole, which is then replaced either by the same dipole, 
or by a pair of dipoles with one common leg. One may naively expect the general Hamil- 
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tonian for dipole evolution to have the same structure as if MSW , but in terms of dressed 
(creation and annihilation) operators. However, even at large N c , the actual dipole evo- 
lution in QCD proceeds in a more subtle way: The bilocal operator p a x p h y within -£/brem 
does not annihilate a dipole, rather it rotates this twice, around different directions in 
color space (see the last line in Eq. (3.15)). These rotations are then compensated by the 
remaining operators in the Hamiltonian (those involving Wilson lines), which can either 
restore back the original dipole, or split it into two new dipoles (cf. Eq. (3.18)). Thus, 
although the initial and final states — prior to, and after the evolution — are colorless 
dipoles, the evolution proceeds through colorful configurations at intermediate steps. 



4 Evolution equation for the dipole densities 

In the previous section, we have expressed the evolution of the dipole picture in two equiv- 
alent ways, which both deal with the ensemble of dipole correlations in the onium: (i) 
as a renormalization group equation for the color glass weight function Z T [p], Eq. (3.22), 
and (ii) as a Master equation for the dipole probability densities P^({zi}; r), Eq. (3.23). 
Alternatively, the same physical evolution can be expressed as an hierarchy of ordinary 
evolution equations for the dipole /c-body densities (k > 1), which in principle can be 
derived from any of the 'functional' evolutions alluded to above provided one knows the 
corresponding expressions for the dipole densities. The derivation based on the Master 
equation (3.23) has been presented in Refs. [4,16]. In what follows, we shall derive the 
same equations from the evolution (3.22) of the color glass. This requires to properly 
identify the operator expressing the dipole number density in the p-representation. 

For more clarity, let us first recall the abstract definition of the dipole number operator, 
where by 'abstract' we mean an operator which is independent of p, and thus of the way 
in which we measure the dipole distribution. The abstract dipole number density operator 
reads (for a iV-dipole configuration) [3, 4] 

N 

n N (x,y; {z,}) = ^^\z z _ 1 - x)S^(z i - y), (4.24) 
and its expectation value n T (x, y) = {n(x, y)) T is obtained as 

oo - 

n T (x,y) = / dT N P N ({z i };T)n N (x,y;{z i }). (4.25) 

N=l J 

By taking a derivative with respect to r in the above equation and using the Master equa- 
tion (3.23), one can check that n T (x,y) obeys the BFKL equation (2.24). Furthermore, 
the abstract operator expressing the dipole pair density is defined as [4] 

N 

n%\x l ,y l ;x 2 ,y 2 ) = Y. * (2) (*i-i - *i)<* (2) (^ - 2/i)<* (2) (^-i - x 2 )6< 2 \z k - y 2 ), (4.26) 

j,k=i 
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where the sum is restricted to different pairs of dipoles since we do not want to count 
the same dipole twice. The equation obeyed by the corresponding expectation value 
nj®(xi, yi, x 2 , y 2 ) has been obtained from the Master equation (3.23) in Refs. [4,16] 
(see Eq. (5.15) in Ref. [4]), and will be rederived below from Eq. (3.22). Abstract /c-body 
dipole density operators n$ with k > 1 can be similarly defined [4, 16]. 

Let us now turn to the representation of these dipole operators in the color glass formalism, 
where a dipole is described as a source for small-rr gluons. Within the Dipole Model, 
where a dipole can radiate only two gluons, the dipole number densities are built with 
the 'bare dipole' annihilation operator introduced in Eq. (3.5). Specifically, by using the 
commutation relation (3.6), one finds (for large N c and x 7^ y) : 

00 N 1 

D (x,y)Z™[p] » £ /dr^P^S^^i-^y + Viv^n^-i.^)^]- 

N=l J i=l Z j^i 

(4.27) 

After averaging over p, which is tantamount to replacing Dq — > 1 (since the functional 
derivatives within D\ give zero after integration by parts), this yields: 

(D (x,y)) T = jD[p}D (x,y)Z™[p] 

00 „ ^ l 

~ dV N Pn{t) 9 (tzi-iv&ziv + S Zi _ iy 5 ZiX ) 

N=l J i=l Z 

= (1/2) (n T (x,y) + n T (y, x)), (4.28) 

where in writing the last line we have recognized the average dipole number density 
according to Eq. (4.24). Note that the measure of the dipole density provided by the 
charge operator D is symmetrized between the quark and antiquark legs of the dipole. 
One can similarly check that, at large N c , 

{D (x 1 ,y 1 ) D (x 2 ,y 2 )} T = - (n^ ) (x 1 ,y 1 ;x2,y 2 ) +n ( ? ) (y 1 ,x 1 ;x 2 ,y 2 ) 

+ n^(x 1 ,y 1 ;y 2 ,x 2 )+n^(y 1 ,x 1 ;y 2 ,x 2 )), (4.29) 

so long as we restrict ourselves to different dipoles; that is, we exclude configurations 
such that the two measured dipoles are identical with each other: {x-± = x 2 ; y x = y 2 } 
or {cci = y 2 ; y\ = x 2 }. Besides, we exclude, as usual, the zero-size dipoles; that is, we 
assume X\ 7^ y\ and x 2 7^ y 2 . 

We thus see that, within the DM, the (bare) dipole annihilation operator plays also the 
role of a number operator. This degeneracy is possible because our present use of creation 
and annihilation operators is somewhat different from their standard use in quantum 
mechanics: Whereas the weight function (3.1) plays naturally the role of a "ground state" 
for the (bare) dipole "Fock space" , the color glass expectation value on this 'ground state' 
involves Z™ itself, and not \Z™\ 2 . 
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The above identification of D with the dipole number operator is also consistent with 
the color glass evolution of the DM, as encoded in Eq. (3.7) for Z^ M . Using the latter 
or, more directly, acting with the MSW Hamiltonian (3.4) on the relevant operators built 
with D , one can indeed check that the correlations of D introduced in Eqs. (4.28) and 
(4.29) obey the evolution equations expected for the respective dipole densities [4,16]. 
The most efficient way to perform this calculation is to use the representation of Hmsw 
as a fourth-order differential operator (cf. Eq. (3.2)) acting on functionals of p a (x). 

Let us now turn to the most interesting case, in which a dipole is allowed to radiate 
arbitrarily many gluons and the respective creation operator involves Wilson lines, cf. 
Eq. (3.8). In the previous manipulations leading to the dipole picture, we did not need 
to introduce the corresponding annihilation operator, because the natural candidate in 
that respect — the bilocal operator p a (x)p a (y) — does not enter the Bremsstrahlung 
Hamiltonian (unlike what happens in the Dipole Model). Yet, as we shall now explain, 
this operator plays an important role also in the general case, not as an annihilation 
operator, but rather as a dipole number operator. Specifically, if one defines 

D(x,y) = -J—p^p^y) for x + y, (4.30) 

9 ■'"c 

which differs from the corresponding DM operator, Eq. (3.5), only by the presence of the 
time argument x + = oo (that we have temporarily reintroduced for more clarity), then 
this operator is gauge invariant and obeys the following commutation relation at large N c 



[D(x,y),D\u,v)] a ^(S ux 5 vy + 5 uy 5 vx ) D\u,v), (4.31) 

which qualifies it as a number operator, as anticipated. To verify the above commutator, 
note that, for x ^ y, Eq. (3.15) yields 

D xy ,Dl v ] = —L(6 W - S xv )tr(t a W u Wl) p a y 
-^- 2 (ty U -tyvMt a W u Wl)p a x 

+ ^ (SuxSvy + SuySvx) AL; (4.32) 

where only the last term in the r.h.s. survives when acting on the onium weight function 
(3.10) and at large N c . (This can be checked through manipulations similar to those in 
Eqs. (3.19)-(3.21).) 

Eq. (4.31) implies that — the weight function for a given configuration of N dipoles; 
see Eq. (3.13) — is an eigenstate of D at large iV c (compare to Eq. (4.27)) : 

1 N 

D(x, y) Z N ({ Ui , v^) « - (dzi^Jziy + S Zi _ iy 5 ZiX \ Z N ({u h (4.33) 

1 i=i 
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which in turn implies (cf. Eq. (4.25)) : 

(D(x, y)) T = J D[p]D(x,y)Z T [p] » \{n T {x, y) + n T (y, xj). (4.34) 

More generally, when acting on Z T [p] with a string of k D-operators corresponding to 
non-identical dipoles, then one measures the /c-body dipole density (symmetrized under 
the exchange of the quark and antiquark legs of each dipole); e.g., for k = 2, one finds 
the analog of Eq. (4.28). Note that, although two operators like D( Xl , t/i) and D(x 2 , 2/2) 
do not commute with each other, this non-commutativity is irrelevant when computing 
their action on the onium weight function and for large N c . 

In the remaining part of this section, we shall verify that the above interpretation of 
the operator D(x,y), Eq. (4.30), as the dipole number operator is consistent with the 
evolution of the correlations of p generated by the Bremsstrahlung Hamiltonian (2.7) 
at large N c . In Sect. 2, we have already performed this check for the 2-point function: 
when acting on p a (x)p a (y), -£/brem generates the BFKL equation (2.24) for (D(x, y)) T , in 
agreement with Eq. (4.34). In what follows, we shall use the results of Sect. 3 to perform 
the corresponding check for the dipole pair density. That is, we shall verify that the 
evolution equation satisfied at large N c by the following 4-point correlation function: 

(£>(*!, yi)D(x 2 , y 2 )) T = -J—( p -( Xl ) p -( yi ) p \ X2 ) p \y 2 )) T: (4.35) 

9 * *c 

is consistent with the known evolution equation for n^(xi,yi, #2,2/2), including the in- 
teresting, 'fluctuation', term (i.e., the term linear in n which describes the formation of 
the dipole pair {(xi,yi), (#2,2/2)} through the splitting of one dipole in the last step of 
the evolution) [4, 16]. 

Using Eq. (3.22), we obtain 

d a 00 r N r r 

^-{D( Xl , yi )D(x 2 , y 2 )) T = -1 ]T dT N P N {r) £ / M m / D[p] D(x u yi)D(x 2 , y 2 ) 
dr 2vr N=1 J i=1 J z J 

x [rf(ui, z) D\z, Vi ) - rtfa, Vi )\ J] D\ Uj , Vj ) 6[p] . 

(4.36) 

To evaluate the r.h.s. of this equation, we need to move the operators DD all the way 
to the right using the commutator Eq. (4.31) and then set = 1. Clearly, the following 
commutator vanishes 

[D xiyi , [D XW2 , Dl kVk \] (any k), (4.37) 

unless the two external dipoles completely overlap, an uninteresting situation that we 
exclude. Moreover, terms like 

[D xiyi ,Di. v .}[D X2y2 ,Dl kVk ] (j,k^i), (4.38) 
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do not contribute either, because the corresponding coefficient vanishes after setting = 
1. Therefore, we can replace the operator part in Eq. (4.36) with 

Dxiy 1 D x . 2 y 2 (Dl l . z Dl, Vi — D^^.) Y\ ^lijVj 
- ~_D xiyi , {Dl iZ Dl Vi - Dl m )] (£ [D X2y2 , Dl jVj ]) + {1 - 2} 

+[D X1V1 , Dl iZ ] [ D X2V2 , D\ Vi \ + {1 <- 2}. (4.39) 

The second line of Eq. (4.39) describes the BFKL evolution of one of the two dipoles in 
the pair (and for large N c ). To see this, note that, e.g., 



D xiyi , (JJ^ VbiZ D\ Vi 



Dt=l 



— — E |-^a;i«iyi^a;itti + M- Ui y xXx & yxVi — ^ M. x\y\z& x^ufi yWi + {^1 
* z 

2 ^ ' J ^-^-xizyi^xiUi^zVi ~l" ■M-zyixi$zu i $yiv i ■M-x\y\z$x\u i $y\v i ~\~ { x l l/l} j" i 



(4.40) 



which represents the BFKL evolution of the dipole cf. Eq. (2.24). 



The last line of Eq. (4.39) is the fluctuation term that we are primarily interested in. 
For large N c , the commutators there can be evaluated according to Eq. (4.31). Note 
that the use of Eq. (4.31) (instead of the exact relation (4.32)) automatically avoids 
contributions in which the external dipoles get mixed with each other under the action of 
the commutators (that is, contributions where the two factors of p from a same external 
dipole get contracted with factors of W coming from different internal dipoles); this was 
to be expected, since such contributions are indeed suppressed at large N c . 

After also integrating over z and summing over N, we finally obtain 



Q T (D Xl ,y 1 D X2i y 2 } T — -f^BFKL + -^BFKL (D Xl>yi D X2) y 2 } 1 



+ 



BFKL 



47T I 



l xiy 2 x 2 



fix 2 yi {-Dxiy2)r 4" •M XlX2 y 1 5y 1 y 2 (D Xl x 2 ) 1 



-\~ ■AAy 1 y 2Xl 5 XlX2 ^\Dy x y 2 ^j T + y\X 2 X\^> X\y 2 



(Dy lX2 ) T j") 



(4.41) 



with the compact notation 

HRFKL{Dx 1 ,y 1 D X2 ^y 2 ) T = —j — M. {x\, y±, z) (D Xl :V1 D X2 :V2 ) T 

+ M(x 1 ,z,y 1 ) {D XuZ D X2 ,y 2 ) T + M(z,y u xi) (D z , yi D X2jy2 ) T . 

(4.42) 

In identifying the expectation values in the r.h.s.'s of the above equations, we have used 
the (dressed-dipole version of the) relations (4.28) and (4.29). The fluctuation terms are 
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the terms proportional to the average dipole number density (D) T in Eq. (4.41). Using 
Eq. (4.28) and (4.29) once again, one can check that Eq. (4.41) is equivalent to Eq. (5.15) 
of Ref. [4]. This confirms the interpretation of {DD) T as the average dipole pair density. 



5 Bremsstrahlung effects on dipole-onium scattering 

So far, we have almost exclusively focused on the evolution of the target wavefunction : 
the large--/V c evolution in the dilute regime has been formulated either as an evolution 
equation for the onium weight function, Eq. (3.22), or as a set of coupled evolution equa- 
tions for those correlations of p which have the meaning of dipole densities (see, e.g., 
Eqs. (2.24) and (4.41)). But the scattering amplitudes for external projectiles involve also 
different correlations of p (or, more precisely, of a = a^(x); see, e.g., Eq. (2.5)), so it 
would be interesting to establish the corresponding evolution equations as well, at least 
at large N c . As anticipated in the Introduction, the evolution equations for scattering am- 
plitudes generated by -£^brem are a priori more general than those previously derived in 
the literature from the dipole picture [4-6] , in the sense of including additional processes 
in which individual target dipoles exchange more than two gluons with the projectile. 
However, our main purpose in this section is not to derive such more general equations, 
but rather to give an argument, based on an explicit example, that the additional effects 
are in fact suppressed at sufficiently high energy. 

The simplest scattering problem which is sensitive to Bremsstrahlung (i.e., to gluon num- 
ber fluctuations) in the target wavefunction is the scattering with two external dipoles. 
(For a single external dipole, -£/brem generates the standard BFKL equation; see below.) 
To lowest order in perturbation theory, each external dipole can exchange two gluons with 
the target, cf. Eq. (2.1). Therefore, the leading-order contribution to the amplitude for 
the simultaneous scattering of the incoming dipoles reads 



(T (2) (^i, Vl ; x 2 , y 2 )) T ~ ((ao(aJi) - a a {Vi)) 2 (a b (x 2 ) - a b (y 2 )f) ^ , (5.1) 



and involves a total exchange of four gluons (two with each external dipole). These four 
gluons can be absorbed either by two different dipoles in the target wavefunction, or by 
a single such a dipole, and the two type of processes are parametrically of the same order 
in a s and 1/N C . In particular, they both contribute in the large- N c limit. Still, as we 
shall argue below, the corresponding contributions behave differently when increasing the 
energy, in such a way that multiple exchanges with the same target dipole are relatively 
suppressed at high energy. 

To see this, we need to evaluate the average in Eq. (5.1) with the onium weight function 
(3.10). The scattering operator for a single (external) dipole, that is, 





(5.2) 
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can be conveniently expressed in terms of the (target) dipole number operator introduced 
in Sect. 4. Namely, after inverting the Poisson equation (2.4) to relate a a to p a , one obtains 

a a (x) - a a (y) = f G(u\x, y) p» , G(u\x, y) = -L In (ti ~ y) * , (5.3) 
J Att (u — x) z 

u 

which then allows us to successively write 



T (x,y) 



4N r 



g(u\x, y)G(v\x, y) p a {u)p a (v) 



8N r . 



G(u\x,y)-G(v\x,y) p a {u)p a (v) 



= J A (x,y\u,v)D(u,v), 



(5.4) 



where in going from the first to the second line we have used the fact that the system is 
globally color neutral : J u p a (u) = 0. The operator D(u, v) has been defined in Eq. (4.30), 
and .4.0 is the amplitude for dipole-dipole scattering in the two-gluon exchange approxi- 
mation and for large iV c : 



A (x,y\u,v) = 



In 



(x — v) 2 (y — u) 2 
(x — u) 2 (y — v) 2 



1 2 



(5.5) 



It is now straightforward to compute the average scattering amplitude for a single exter- 
nal dipole: using Eq. (4.33) together with the symmetry of the dipole-dipole amplitude 
Aq(x, y\u, v) under the exchange u <-> v, one immediately obtains 



(T(x,y)) T = J Ao(x,y\u,v)n T (u,v), (5.6) 

which is the relation expected within the dipole picture [1-3]. By using this relation 
together with Eq. (2.24) for n T , one can show [3] that (T(x,y)) T obeys the standard 
BFKL equation, as anticipated. 

To similarly evaluate Eq. (5.1), we shall use Eq. (5.4) for both external dipoles. We have 



(T ( - 2 \x 1 ,y 1 ;x2,y 2 ))T = J Ao(x 1: yi\u 1: «i) Ao(x 2 , y 2 \u 2 , v 2 ) (D(ui, v{) D(u 2 , v 2 )) T , 

Ui,Vi 

(5.7) 



where after using again Eq. (4.33) and the symmetries of the tree-level amplitude ^4o, 
one can effectively replace (at large iV c ) 
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N N 

(D(u u Vl ) D(u 2 , v 2 )) T -> ( £ 5 (2) (^-i - u^ 2 \z t - Vl ) £ ^ (2) (^-i " u 2 )5^ 2 \z k - v 2 ) 

i=l k=l 

= n {2 \u u vi, u 2 ,v 2 ) + 5 (2) (u 1 - u 2 )5 (2) (v t - v 2 )n T (u 1: vi). 

(5 

We have recognized here the average dipole density in the target n T and also the average 
dipole pair density according to their definitions in Eqs. (4.24)-(4.26). When the 
above expression is inserted in the r.h.s. of Eq. (5.7), we finally obtain: 

(T (2) (a3!, yi, x 2 , y 2 )) T = J A (x 1 , yi|«i, v ± ) A {x 2 , y 2 \u 2 , v 2 ) n[ 2) («i, v t ; u 2 , v 2 ) 

Ui,Vi 

+ J A (x 1 ,y 1 \u,v)A (x 2 ,y 2 \u,v)n T (u,v), (5.9) 

with a clear physical interpretation for the two terms in the r.h.s. : The first term, propor- 
tional to rS 2 \ describes the scattering with two different dipoles in the target, while the 
second term, proportional to n T , represents a four-gluon exchange with a single dipole. In- 
cidentally, the above derivation of Eq. (5.9) confirms that no ordering ambiguities appear 
in the evaluation of the scattering amplitudes at large iV c . 

Previous applications of the dipole picture to scattering [1-4, 6] were all based on the 
assumption that a dipole can exchange only two gluons, and thus they have ignored the 
second term in Eq. (5.9). Using only the first term there, together with the known equation 
for n[ 2 ) (essentially, Eq. (4.41)), one has derived an evolution equation for (T^) T [4,6] 
which includes the essential 'fluctuation term' through which (T^) T gets built from (T) T 
in the dilute regime. (This is induced by the fluctuation terms in Eq. (4.41).) In principle, 
one can similarly use Eq. (5.9) together with the known equations for and n T to 
deduce the most general evolution equation for {T^) T in the dilute regime and at large 
N c . In practice, this might be however tedious, as it requires to invert Eqs. (5.6) and (5.9) 
in order to express and n T in terms of the scattering amplitudes. 

But even without any detailed calculation, it is clear by inspection of Eq. (5.9) that at 
high energy (while still in the dilute regime, though, for the dipole picture to apply), the 
contribution involving the dipole pair density rS 2 ^ is in fact the dominant one: Indeed, 
with increasing energy, n T grows like a BFKL pomeron, n T ~ exp{tua s r}, whereas 
grows like a Pomeron squared: ~ exp{2cua s r} (tu is a pure number). Therefore, the 
contribution proportional to dominates as soon as a s r ^ 1, which is the interesting 
regime at high energy 5 . This implies not only that, for a s r ^ 1, one can neglect n T next 
to in the r.h.s. of Eq. (5.9), but also that, in writing an equation for (T^) T which 
should describe its evolution from the low energy regime at r ~ up to the high energy 
where a s r ^ 1, it is enough to keep trace of the terms coming from the evolution of . 
This is what has been done in Ref. [4,6]. 

5 This condition leaves a parametrically large window for the applicability of the dipole picture, 
which ceases to be valid when n T ~ l/a 2 s or a s r ~ ln(l/a;?) [1,4]. 
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Although certainly correct (in view of the previous arguments), the last statement is 
nevertheless quite subtle, as it can be appreciated when trying to identify the fluctuation 
term (the contribution proportional to (T) r ) in the evolution equation for (T^) T derived 
from Eq. (5.9). By taking a derivative w.r.t. r in Eq. (5.9) and using the equations for 
and n T , one finds two type of contributions which are proportional to n T (and thus to (T) T , 
cf. Eq. (5.6)) : the 'genuine' fluctuation term coming from the evolution (4.41) of nf\ 
and the terms describing the BFKL evolution (2.24) of n T . Both types of contributions 
are parametrically of the same order — namely, of order a s a 2 s {T) T — , and thus contribute 
on equal footing to the growth of (T^ 2 ^) T in the dilute regime. This seems to contradict 
the results in Ref. [4,6], where only the 'genuine' fluctuation term has been included. 

However, this contradiction is only illusory: First, to properly identify the fluctuation 
terms, one need to reexpress everywhere in terms of (T^) T , according to Eq. (5.9); this 
operation modifies the fluctuation terms since, schematically, i^nf = {T^) T — Aln T . 
Second, the additional fluctuation terms (besides the 'genuine' one) which persists after 
the previous operation have only the role to compensate the spurious double-Pomeron 
contribution which emerges from the evolution of that piece of the initial condition (T^) 
which is introduced by the second term in Eq. (5.9). Therefore, if one starts with the 
following initial conditions at r = (i.e., a target made with a single dipole): 

n (u,v) = 5 {2) (u-u )5 {2) (v -v ), nf ] = , (5.10) 

which in turn implies 

(T (2) (a;i, yi, x 2 , y 2 ))o = Ao(x u I/i|«o, «o) Ao(x 2 , 2/ 2 |«o, v Q ) , (5.11) 

then the dominant contribution to (T^) T at high energy — the one which grows like a 
double-Pomeron, ~ exp{2ua s T} — comes entirely from the 'genuine' fluctuation term, 
and not from the evolution of the above initial condition. The latter will give only a 
subleading contribution, ~ exp{ud s r}. 

These considerations can be clarified with the help of a simple example, borrowed from 
Ref. [4], which has the advantage to be easily solvable while keeping the non-trivial 
features of interest here 6 . Namely, let us replace the dipoles par point-like particles 
which live at a fixed point (so there is no spatial dimension involved in their dynamics), 
and whose number distribution evolves according to the following hierarchy of equations: 



dn T dn^ 



= an T , — ^— = 2a 



(5.12) 



dr dr 

which mimic the equations satisfied by the dipole densities in the dipole picture. The 



6 It is likely that the following argument is similar to that developed in a different context by 
Braun and Vacca [24], but we have not been able to clearly establish the correspondence between 
the two problems. See also the discussion in Ref. [6]. 
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solution corresponding to the initial conditions n(0) = n and n( 2 )(0) = reads: 



n T = n exp(ar), 

= 2no exp(2ar) — 2no exp(ar), . . . (5.13) 

where has been generated by the term linear in n T in the r.h.s. of the second equation 
(5.12); this is the analog of the 'genuine fluctuation term' in the present model, and plays 
the role of a source for (so like the actual fluctuation terms in Eq. (4.41)). If one 
further introduces the analog of the 'scattering amplitudes' : 

T T = n T , T T (2) = An 1 *? + Bn T , (5.14) 

then clearly 

T^ 2) = 2An [exp(2ar) - exp(ar)] + Bn exp(ar) w 2An exp(2ar), (5.15) 

where the approximate equality holds at large time (ar 1), and the corresponding 
contribution to T^> is entirely coming from (like in Eq. (5.9)). Let us now construct 
the evolution equation for and follow the fluctuation terms: 

— — = A—^— + B— - = 2Aa n?> + n T + Ban T 
dr dr dr 1 J 

= 2aT r (2) + (2A-B)aT T , (5.16) 

where the second line identifies (2A—B)a T T as the fluctuation term. This involves the two 
types of contributions alluded to before: 2 Act T T is the genuine fluctuation term introduced 
by the evolution of whereas (—B)aT T has been induced by the BFKL evolution of 
n T . With T T = n exp(o;r), the above equation is solved by 

T (2) =T (2) exp ( 2ar ) + (2A - B)n [exp(2ar) - exp(ar)], (5.17) 

where both types of fluctuations seem to contribute to the dominant, 'double-Pomeron', 
rise at large time. However, by recalling that = Bn , cf. Eq. (5.14), one immediately 
sees that the 'double-Pomeron' terms proportional to B do actually cancel between the 
contribution of the initial condition and that of the fluctuation terms. The remaining 
'double-Pomeron' term, proportional to A, is the one generated by the genuine fluctuation 
term, so like in the direct calculation leading to Eq. (5.15). Thus, the same dominant 
behaviour at large time would have been obtained by solving the simplified equation: 

dT^ 

= 2aT^ + 2AaT T , (5.18) 

dr 

(which arises by assuming that = An^) together with the initial condition T^ 1 = 0. 
Eq. (5.18) is the analog of the equation for (T^) T derived in Refs. [4, 6], whereas Eq. (5.16) 
corresponds to the more complete equation that would be obtained in QCD at large N c 
after including the effects of Bremsstrahlung in the target wavefunction. 
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The simple example above also emphasizes the importance of correctly adjusting the 
initial condition to the approximations that we perform on the evolution equation: If, 
as in Refs. [4,6], we restrict ourselves to the two-gluon exchange approximation in the 
construction of the evolution equation for (T^) T , then the same approximation must be 
performed also on the initial condition. For instance, if at r = the target is a bare 
dipole, cf. Eq. (5.10), then the initial condition to be used for the equations in Refs. [6] 
is (T^)q = 0, and not Eq. (5.11). This amounts to consistently neglect contributions like 
the second term in Eq. (5.9) at all places. 



6 Conclusions 

In this paper we have shown that the recently developed theory for Bremsstrahlung in the 
QCD evolution with increasing energy [8-10] is consistent with the dipole picture in the 
large- A^ c limit. The characteristic feature of this theory is that the quark and the antiquark 
parts of a color dipole are represented as color sources which can radiate arbitrarily many 
gluons with relatively small longitudinal momenta. This generalizes previous 'color glass' 
descriptions of the onium wavefunction, in which the individual dipoles were allowed to 
radiate only two gluons [3-6]. The energy evolution of the ensemble of color sources is 
known for arbitrary N c — this is described by a two-dimensional, Hamiltonian, field 
theory with SU(iV c )-like Poisson brackets — , and our objective in this paper has been 
to demonstrate that, for large N c , this evolution can be reformulated in terms of color 
dipoles which evolve through dipole splitting. To that aim, we have identified the operator 
which describes, within the Hamiltonian theory, a quark-ant iquark dipole 'dressed' by the 
radiation, and then we have shown that, at large N c , the action of the Hamiltonian on a 
collection of such dipoles consists in the splitting of any of the original dipoles into two new 
dipoles with one common leg. This confirms the fact that the effective degrees of freedom 
for high energy evolution in QCD at large N c and in the dilute regime are quark-antiquark 
color dipoles, as originally demonstrated at an abstract level (i.e., without specifying the 
way how the dipoles are actually measured) in the pioneering papers by Mueller [1, 2]. 

Whereas the emergence of the dipole picture at large N c was to be expected (in view of 
the general results in Ref. [1]), the technical manipulations necessary to demonstrate it 
turned out to be quite complex. In particular, we have found that the structure of the 
Hamiltonian cannot be further simplified when going to the large-iV c limit. Rather, it is 
its action on the 'onium weight function' (built with specific dipole creation operators) 
which reduces to dipole splitting at large N c . Thus the present construction provides 
a rather subtle generalization of previous formalisms based on the two-gluon exchange 
approximation [3, 5]: Whereas the structure of the wavefunction in terms of dipole creation 
operators is formally the same (with different meanings for the creation operators though), 
the corresponding Hamiltonians are very different, and cannot be directly related to each 
other via some large- N c approximations. 

This situation is reminiscent of that encountered in a recent large-iV c analysis of the 
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JIMWLK evolution in the high density regime [7], and in fact the duality between the 
two Hamiltonian field theories (JIMWLK and Bremsstrahlung) [9, 10] has played an im- 
portant role in the present analysis. Namely, the action of the JIMWLK Hamiltonian on 
scattering operators for external dipoles turns out to be dual to that of the Bremsstrahlung 
Hamiltonian on creation operators for internal dipoles. Because of that, many of the tech- 
nical manipulations in this paper are similar (more properly, dual) to those encountered 
in the derivation of the Balitsky equations [22] from the JIMWLK equation [13,20]. 

We have finally investigated the effects of Bremsstrahlung on the scattering between two 
external dipoles and a dilute onium. We have shown that, although some new effects 
appear, as associated with the double scattering off a same target dipole, such effects are 
in fact suppressed at high energy as compared to those already included in the two-gluon 
exchange approximation. This conclusion, which is corroborated by the recent analysis in 
Ref. [18], implies that the evolution equations with 'pomeron loops' derived in Refs. [4-6] 
are indeed the correct equations in QCD at large N c and for sufficiently high energy. 

It remains as an interesting open problem to consider the generalization of some of the 
results obtained here (in particular, in relation with the scattering problem) to arbitrary 
N c . Although the corresponding evolution Hamiltonian is known, and so is also its action 
in terms of Poisson brackets (cf. Sect 2), we expect the corresponding analysis to be 
complicated by the issue of the ordering of the operators. This is a new type of problem, 
which had not been encountered before in the framework of the color glass formalism, 
and so far it is not even clear whether this formalism can be extended to account for 
the non-commutativity of the color charges (see, e.g., the discussions in Refs. [8-11, 18]). 
One may however expect to be able to establish more explicit connections to previous 
approaches within perturbative QCD, which are aiming at the direct calculation of the 
2 — > n gluon vertices at high energy and for arbitrary N c [24-27]. 
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